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A Young White Dwarf with an Infrared Excess 
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ABSTRACT 

Using observations of Spitzer/IRAC, we report the serendipitous discovery of 
excess infrared emission from a single white dwarf PG 0010+280. At a tempera¬ 
ture of 27,220 K and a cooling age of 16 Myr, it is the hottest and youngest white 
dwarf to display an excess at 3-8 /irn. The infrared excess can be fit by either an 
opaque dust disk within the tidal radius of the white dwarf or a 1300 K black- 
body, possibly from an irradiated substcllar object or a re-heated giant planet. 

PG 0010+280 has two unique properties that are different from white dwarfs 
with a dust disk: (i) relatively low emission at 8 /irn and (ii) non-detection of 
heavy elements in its atmosphere from high-resolution spectroscopic observations 
with Keck/HIRES. The origin of the infrared excess remains unclear. 

Subject headings: circumstellar matter - minor planets, asteroids: general - white 
dwarfs 


1. Introduction 


Radial velocity and transit surveys show that extrasolar planetary systems are prevalent 
around main-sequence stars in the Milky Way. Such direct evidence for planets orbiting white 
dwarfs is sparse although a ~ 7Mj mass object has been detected with direct imaging at 
2500 AU from the white dwarf WD 0806-661 (Luhman et al. 2011). There exist strong 


indirect arguments that planetary systems are common around white dwarfs (Zuckerman 


et al. 

2003; 

Koester et al. 

2014; 

Jura & Young 

2014) 


Theoretical calculations show that 
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planets and minor planets beyond a few AUs can survive the red giant stage of the star 
and persist into the white dwarf phase (Burleigh et al. 2002 Jura 2008). The orbits of 


these planets expand correspondingly and this dynamical rearrangement causes additional 
instability in the system. Some surviving minor planets in the system get perturbed, enter 
into the tidal radius of the white dwarf, and are eventually disrupted and accrete onto the 


white dwarf (e.g., 

Debes & Sigurdsson 2002 

Debes et al. 2012 

. These tidal disruption events 

can continue over a few Gyr of the white dwarf cooling (e.g., 

Bonsor et al. 

2011 

Veras et al. 


2013) and produce two observables: (i) a disk within the tidal radius - both the dust and 
gas component have been identified ( Gansicke et al.||2006 Xu &; ,Jura||2012 ): (ii) atmospheric 
enrichment of heavy elements in addition to a white dwarf’s primordial hydrogen or helium 
atmosphere - 25-50% white dwarfs show heavy elements from accretion of circumstellar 


material (Zuckerman et al. 2003; Koester et al. 2014) 


The dust disk occurrence is about 4% for white dwarfs between 9500 K and 25,000 K 
(Barber et al. 2012 Rocchetto et al.|2014 ). Until now, all dusty white dwarfs also have heavy- 
element-enriched atmosphere, which provides unique information about the composition of 
its accreting material (e.g., Klein et al. 2010[ Gansicke et al. 2012 Xu et al. 2014). In 


addition, there have been a lot of efforts to look for low-mass and substcllar companions 
around white dwarfs by searching for infrared excess (e.g., Zuckerman & Becklin 1987). 


There are a handful of white dwarf - brown dwarf pairs and the frequency is less than 0, 


(Farihi & Christopher 2004 Farihi et al. 2005). The main difference between a dust disk 


and a companion is the color of the infrared excess. 


In this Letter, we report the serendipitous discovery of infrared excess around PG 0010+280 
- the hottest and youngest hydrogen-atmosphere white dwarf with a 3-8 /mi excess. We also 
put a strong constraint on the mass accretion rate. 


2. Observations and Data Reduction 


PG0010+280 was discovered as a hydrogen-atmosphere white dwarf in the early 80s 


(Wegner 1983). Recently, Gianninas et al. (2011) analyzed this star and derived T* = 


27,220 K, log g = 7.87 and d= 152 pc. This corresponds to a white dwarf mass of 0.57 M 0 
and main sequence mass of 1.8 M 0 (Williams et al. 2009]). The main sequence lifetime is 
about 2.3 Gyr and the white dwarf cooling age is 16 Myi[]] (Holberg & Bergeron 2006). 


^rom the cooling models presented in 


http: //www. astro.umontreal.ca/~bergeron/CoolingModels 
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2.1. Spitzer/IRAC and WISE 


PG 0010+280 was observed with all four IRAC bands in 2006 in program 30856. The 
observation was performed by using a 30 sec frame time in the setting of a random nine-point 
large-amplitude dithering. The total on target time was 270 sec in each IRAC channel. 


Following data reduction procedures outlined in Jura et al. (2007a) and Xu & Jura 


(2012), we used MOPEX to process each exposure and combine them into a final mosaic 
with a subsample of 2 in native pixels (i.e., 076 pixel _1 ). PG 0010+280 has a proper motion 


of 07354 in RA and -07116 in DEC according to the PPMXL catalog (Roeser et al. 2010) 
and it is located at 00:13:21.09 +28:20:19.70 in the IRAC data. As shown in Figure [lj 
PG 0010+280 is in a clean field and well detected in the first three IRAC bands. IRAF was 
used to perform aperture photometry with the following parameters: (i) an aperture radius 
of 2 pixels (172); (ii) a sky annulus of 12 - 20 native pixels (1474 - 2470). We applied aperture 
correction according to the IRAC instrument handbook. There are two parts contributing 
to the uncertainty, (i) 5% calibration uncertainty (Farihi et al. 2008a) and (ii) measurement 
error. The final values are 145 ± 8 /rJy, 88 + 6 /iJy and 89 ± 21 /iJy for IRAC-1, IRAC-2 
and IRAC-3 band, respectively. In IRAC-4, PG 0010+280 is only marginally detected. We 


followed procedures described in Farihi et al. (2008a) and find a 3cr upper limit of 36 /rJy. 


PG 0010+280 was also observed with WISE between 2010 and 2011. In ALLWISE 


data release (Cutri et al. 2013), it is well detected in WISE bands 1 and 2 with a flux of 


136 ± 6 /rJy and 92 ± 10 /uJy, respectively, comparable with the IRAC fluxes. The WISE 
position at 00:13:21.10 +28:20:19.60 is consistent with the IRAC position after correcting 
for its proper motion. This shows that the infrared excess is physically associated with 
PG 0010+280 rather than from a background source. 


2.2. UKIRT 


In 2014, PG 0010+280 was observed with the Wide Field Camera (WFCAM) on the 
United Kingdom Infrared Telescope (UKIRT) in the J, H and K bands. We used a 5 point 
jitter with a frame time of 5 sec for J band and 10 sec for H and K bands. A few exposures 
were taken in each band and the total on target time was 75 sec, 150 sec and 250 sec for J, 
H and K band, respectively. The data were processed by the pipeline from the Cambridge 
Astronomical Survey Unit (CASU). The final images are shown in Figure [I] CASU pipeline 
also performed flux calibration following the procedure developed for the UKIRT Infrared 


Deep Sky Survey (UKIDSS) (Hewett et al. 2006 Hodgkin et al. 2009). The fluxes in indi¬ 


vidual exposures are consistent with each other and the averaged values are 460 ± 4 /iJy, 
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300 ± 3 /iJy and 191 ± 3 pJy for J, H and K band, respectively. In all three bands, the 
star image looks symmetric and the FWHM is about 2 pixels, corresponding to 078. At a 
distance of 152 pc, the infrared excess should come within 60 AU of PG 0010+280. 


2.3. Keck/HIRES 


High-resolution spectroscopic observations of PG 0010+280 were performed with the 
HIRES (Vogt et al. 1994) on the Keck I telescope. The red and blue collimator were used 


separately on the nights of September 17th and 20th, 2013 under good weather conditions. 
The total on target time was 4800 sec and 1800 sec in the red and blue, respectively. On 
both nights, the C5 slit was used with a width of T/148, making the resolution around 40,000. 
The flux standard Feige 110 was observed on both nights. 


Following Klein et al. (2010) and Xu et al. (2014), we used the MAKEE software to 


locate the trace and extract the spectra. Then we used routines in IRAF to correct for 
atmospheric extinction and derived the sensitivity function for each echellc order using the 
flux standard Feige 110. Because the goal of the observation is to detect narrow white dwarf 
photospheric features, the uncertainty in the absolute flux calibration is not relevant. The 
final spectrum has a coverage from 3100 to 9000 A with a few small gaps between echelle 
orders and it looks like a “normal” hydrogen white dwarfs. The spectrum around selective 
elements is shown in Figure [2] 


3. Discussion 
3.1. SED Fits 


We calculated a pure hydrogen NLTE white dwarf model atmosphere using TLUSTY 


& SYNSPEC( 

Hubeny & Lanz 

1995; 

Lanz & Hubeny 

2003 

) with the following parameters: 

T=27,220 K, log g=7.87 ( 

Gianninas et al. 

2011 

). Because 

we are most interested in iden- 


tifying infrared excess associated with the white dwarf, we scale its photospheric flux to 
UKIRT J band. As shown in Figure [3j the observed flux in the infrared is much higher 
than predicted for the white dwarf’s photosphere. We first fit the infrared excess with a 
flat, opaque disk model proposed in Jura (2003). The inner disk radius R ira varies between 
10 R„,rf and 50 R w d and outer disk radius R ou t = 40 R w d to 100 R w d both in steps of 1 R w d 
and the disk inclination cosi increases from 0 to 1 in steps of 0.01. We defined a y 2 with 
5 degrees of freedom, including flux differences between PG 0010+280 and the model at K 
band and IRAC 1-4 bands. 
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There is no unique fit to the SED due to the degeneracy of the disk inclination and 
surface area (e.g., Jura et al. 2007a). We present the model with the smallest y 2 , as shown 
in Figure [3j The required dust temperatures are 1790 K and 960 K at the inner and outer 
radius, respectively. However, the model can not reproduce the high flux around IRAC-1 
and W-l and the low flux at IRAC-4. In other dusty white dwarf systems, the IRAC-4 
flux is typically higher or comparable to the IRAC-3 flux due to the adjacent 10 /am silicate 
emission feature (Jura et al. 2007a, 2009), as shown in Figure |4j It is puzzling that the 
IRAC-4 flux for PG 0010+280 is lower than the flux from a flat disk model. 

The infrared excess around PG 0010+280 can also be fit with a 1300 K blackbody, as 
shown in Figure [3j At a distance of 152 pc, it requires a radius r p of 9 x 10 9 cm, about 1.3 
radius of Jupiter Rj, which is possible for a substellar object. To explore this possibility, 
we fit the infrared excess with different models developed for brown dwarfs and exoplanet 
and find the best fit is the AMES-DUSTY model with T p =1300 K, log g=4.0 and r p = 1.3 tj 


(correspondingly a mass of 10 Mj, Allard et al. 2001). We defer the further discussion to 
section 3.3. 


3.2. Non-detection of Heavy Elements 


No absorption lines from heavy elements were detected in the HIRES spectrum of 
PG 0010+280. To derive the upper limits, we used the same atmospheric model calculated 
from TLUSTY & SYNSPEC (Hubeny & Lanz 1995 Lanz fe Hubeny| 2003) with atomic 


data for individual lines from the Vienna Atomic Line Database (Kupka et al. 1999). The 


computed lines were convolved with a Gaussian function to represent the instrument broad¬ 
ening. The upper limits are listed in Table [l] and some representative spectra are shown in 
Figure [2} At such a high stellar temperature, there is no convection zone and the usual def¬ 
inition for settling times becomes very arbitrary (Koester et al. 2014). Instead, we calculate 
the accretion rate of each element in Table [lj following Gansicke et al. (2012) and assuming 
the system is in a steady state (D. Koester, private communication). 

So far, PG 0010+280 is the only white dwarf with a 3-8 /irn excess but no detection of 
heavy elements. However, it is not entirely unexpected. The optical line strength decreases 
significantly with increasing stellar temperatures (Koester et al. 2014). For PG 0010+280, 
the magnesium upper limit is the most stringent and we assume magnesium is 15% of the 
total mass of the material - the average number in other dusty white dwarfs (Jura & Young 
2014). The upper limit of the total accretion rate would be 1.2 x 10 8 g s _1 , which is 


relatively low given dusty hydrogen-atmosphere white dwarfs have an average accretion rate 
of 6 x 10 8 g s 1 (|Gansicke et al.|2012||Xu et al.|2014). Future observations, particularly with 
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HST/COS which is more sensitive to trace elements, would reveal whether heavy elements 
are present in the atmosphere. 


3.3. Origin of the Infrared Excess 

3.3.1. A dust disk? 

If the infrared excess comes from a dust disk, PG 0010+280, at a temperature of 
27,220 K, would be the hottest dusty white dwarf. The previous record is WD J1537+5151 


at 24,900 K (Barber et ah 2014). Because the disk is opaque, due to mutual shielding, the 


grains reach a much lower temperature compared to a blackbody. The sublimation radius 


can be calculated using equation (1) in Jura (2003). Assuming a sublimation temperature 
of 1200 K, the sublimation radius is 38 R w d, which is well within the tidal radius. The 
innermost region of the dust disk is directly exposured to the starlight and the temperature 


can be much higher. Rafikov & Garmilla (2012) studied this particular region and found 


that a higher dust temperature is possible in the high metal pressure. It is possible to have 
a dust disk at a temperature of 27,220 K for PG 0010+280. 

The presence of dust disks around young white dwarfs is favored by dynamical simula¬ 
tions, which predict a decline of tidal disruption event as the white dwarf cools. For example, 


Debes et ah (2012) found that the peak of accretion happens at a white dwarf cooling age 


of 30 Myr and stellar temperature of 24,000 K. The degree of pollution is also observed to 
decrease as the white dwarfs age and only two dusty white dwarfs have a temperature lower 


than 10,000 K (Bergfors et al. 2014). 


3.3.2. A substellar companion? 

The SED of PG 0010+280 can also be fit with a substellar object at 1300 K, 1.3 r j and 
within 60 AU of the white dwarf. There are a few possibilities. 


A cool brown dwarf. However, a radius of 1.3 ij is too big for a brown dwarf that has 
been cooling for 2.3 Gyr. For example, Allard et al. (2011) derived a radius of 0.86 r j 
for such a brown dwarf, which is about 50% smaller than required to fit the SED. 


A re-heated brown dwarf. A brown dwarf might accrete a substantial amount of mass 
during the red giant stage and would appear to have a similar age as the white dwarf 
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cooling age (Livio & Soker 1983). In this scenario, the substcllar object would be 


3000 K, which is too hot to fit the SED. 

An irradiated substellar object. Assuming the temperature of the object is in thermal 
equilibrium with the stellar radiation, it would be located at 220 R^ (0.015 AU) to 
be 1300 K. The orbital period is 21 h. The substellar object would be outside of the 


tidal radius, which is typically 100 (Jura 2003). A narrow Ha emission core can 


be visible during part of the orbital period, as has been detected around WD 0137-349 


and its companion (Maxted et al. 2006 Burleigh et al. 2006). Ha was covered in our 


HIRES/red setting in two consecutive exposures, each 40 minutes long. No emission 
core was detected. However, we do not have enough time coverage to rule out this 
possibility. 

A giant planet. At an age of 2.3 Gyr, a planet would be too cool to be detectable. 


A re-heated giant planet. Spiegel & Madhusudhan (2012) studied the evolution of 


giant planets around post-main-sequence stars and found they can be significantly re¬ 
heated due to accretion of stellar material. At a radius of 1.3 rj, it is similar to the 


values derived for inflated extrasolar giant planets (Batygin & Stevenson 2010). This 


recent accretion would leave a huge amount of dust in the atmosphere of the substellar 
objects - the prerequisite for AMES-DUSTY model. While for other typical substellar 
objects at 1300 K, the best fit model has little dust contribution (Allard et al. 200l]). 

For the potential substellar object to remain “hot” and display significant infrared 
excess at the white dwarf cooling age t coo i of 16 Myr, it needs to accrete at least the 
mass of AM, 


‘-pi 


A M p = 


tcooi X Lp X Tp 


GM n 


( 1 ) 


where G is gravitational constant. M p , L p and r p represents the mass, luminosity and 
radius of the substcllar object. Assuming its effective temperature T p has remained 
the same since the start of white dwarf cooling, it can be calculated as, 


L P = &SBTp x 47 t = 1.3 x 10 29 erg s 1 (2) 

For a substcllar object as massive as 10 Jupiter mass, it needs to accrete at least 
6 x 10 29 g of material to remain hot after 16 Myr. This corresponds to ~ 3% of its 
own mass and ~ 0.02% of the total mass loss by the 1.8 M 0 main-sequence progenitor 
star. Depending on the location of the object and the mass loss scenario, this hypothesis 
could be possible. Due to possible accretion of carbon-rich material, the composition 



























of the object’s atmosphere could be substantially modified and display CH 4 and CO 
(Spiegel & Madhusudhan 2012). Spectroscopy in the near-infrared might detect some 
signatures from a substcllar object and reveal the nature of the infrared excess. 


4. Conclusions 

We report the discovery of infrared excess around PG 0010+280 - the hottest and 
youngest white dwarf displaying a 3-8 /mi excess. Its unique infrared color as well as the 
non-detection of heavy elements with high-resolution spectroscopic observations suggest a 
possible alternative origin than white dwarfs with infrared excess from a circumstellar disk. 
From fitting the SED, we can not exclude models of either an opaque dust disk within 
the tidal radius or a substcllar object at 1300 K, from an irradiated object or a re-heated 
planet. Future observations, particularly with spectroscopic observations in the ultraviolet 
and near-infrared, could reveal the nature of the infrared excess. 
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partnership among the Caltech, the University of California and NASA. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. The authors 
also wish to recognize and acknowledge the very significant cultural role and reverence that 
the summit of Mauna Kea has always had within the indigenous Hawaiian community. We 
are most fortunate to have the opportunity to conduct observations from this mountain. This 
material is based upon work supported by the National Science Foundation under Grant No. 
1154375. 





Table 1: Derived Abundance Upper Limits for PG 0010+280 


Element 

Line 

EW (mA) 

log n(Z)/n(H) 

Flnx a 
(g s -1 ) 

C 

C 11 4267 b 

< 60 

< -4.5 

< 6.1 x 10 7 

O 

O 1 7772 

< 25 

< -3.5 

< 2.4 x 10 9 

Mg 

Mg 11 4481 6 

< 38 

< -5.5 

< 1.8 x 10 7 

Si 

Si II 6347 

< 22 

< -5.3 

< 2.7 x 10 7 

Ca 

Ca II 3933 

< 40 

< -5.5 

< 4.9 x 10 7 

Fe 

Fe II 5260 

< 14 

< -4.1 

< 1.6 x 10 9 


Upper limit to the accretion rate of each element, based on the formalism described 


Gansicke et al. (2012) (D. Koester, private communication). 


This line is a doublet. 
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Fig. 1.— From top left to bottom right, images for PG 0010+280 in H and K bands from 
the UKIRT as well as four Spitzer/IRAC bands. North is up and East is left. The field of 
view is 30 / . , 0 by 30'/0. Red circles are centered on PG 0010+280 with a radius of 6"0. The 
contrast is adjusted for good viewing and not constant in different images. PG 0010+280 is 
readily detected in H, K and IRAC bands 1-3 and only marginally detected in IRAC-4. 
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Fig. 2.— Keck/HIRES observation of PG 0010+280. The black line represents the data and 
the red line represents the model with a fit to the upper limit (see Table [I]). Wavelengths 
are given in air and in the star’s reference frame. 
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Fig. 3.— Fits to the SED of PG 0010+280, including data from SDSS, UKIRT, 1RAC 
bands 1-4 and WISE bands 1-2. Top: the black dashed line represents the computed white 
dwarf model spectra flux calibrated to J band. Bottom: the infrared excess - the white 
dwarf photospheric flux has been deducted from the data. Three model fits are presented, 
including (i) a dust disk (dashed magenta line) with R, m = 23 R wc i (T* n = 1790 K), R out = 
53 R , w d (T out = 960 K) and cos i = 0.13; (ii) a blackbody (dashed black line) at 1300 K with 
a radius of 1.3 rj; (iii) an AMES-DUSTY model (grey line) for a giant planet (Allard et al. 


2001 ). 
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Fig. 4.— Spitzer/IRAC colors of all white dwarfs with infrared excess with available mea¬ 
surements. PG 0010+280 stands out by having a very small (IRAC-3 - IRAC-4) value. While 
GD 362 has unusually large 10 jim silicate emission feature (Jura et al. 2007b) and therefore 
a big (IRAC-3 - IRAC-4) value; WD 1041+092 and WD 0843+516 also have a high IRAC-4 
flux possibly with a similar origin. 
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